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Abstract. The MAGIC 7-ray observatory has re- 
cently been upgraded by a second 17 m-diameter 
imaging atmospheric Cherenkov telescope at a dis- 
tance of 85 m from the first one. Simultaneous obser- 
vation of air showers with the two MAGIC telescopes 
(stereoscopic mode) will improve the reconstruction 
of the shower axis and solve the ambiguity in the 
impact point occurring in single-telescope mode. 
Also, the stereo observation will result in a better 
angular resolution, energy estimation and cosmic- 
ray background rejection. It is expected that the 
sensitivity of MAGIC improves significantly over 
the full energy range (60GeV - 20TeV). Here, we 
present the performance estimated from Monte Carlo 
simulations. 
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I. The MAGIC Telescopes 

The MAGIC telescopes are a system of two 17 m- 
diameter imaging atmospheric Cherenkov telescopes 
(IACT) used for very high energy (VHE) 7-ray astron- 
omy. The first telescope, MAGIC-I, is in operation since 
2004. The second one, MAGIC-II, was built in 2008 
at a distance of 85 m from MAGIC-I. It is currently 
under commissioning and scientific observation should 
start within a few months Q]. 

The structure of the two telescopes is almost identical, 
based on a reinforced carbon fiber frame supporting 
236 m 2 of mirrors. The two cameras, however, are quite 
different. MAGIC-I has a 3.5° field of view camera 
composed of two types of pixels with different sizes. 397 
small pixels (0.1° diameter) cover a central hexagonal 
region and 180 larger pixels (0.2°) cover a surrounding 
ring. MAGIC-II has an improved camera [2| equipped 
uniformly with 1039 small pixels (0.1°) with a better 
photo-detection efficiency, and covering a field of view 
of 3.5 °. The trigger region is limited to a central region 
of 2 diameter for MAGIC-I and 2.5 for MAGIC- 
II. Both telescopes have a 2 GSample/s digitization data 
acquisition system. 

The two MAGIC telescopes can be operated indepen- 
dently or in stereoscopic mode. The stereoscopic ob- 



servation mode allows a more precise reconstruction of 
the shower parameters as well as a stronger suppression 
of the hadronic showers and other background events. 
Here we present the performance in stereoscopic mode 
(MAGIC phase-II) and compare it to the performance of 
MAGIC in phase-I using only the MAGIC-I telescope. 

II. Monte Carlo Simulations 

The performance of the MAGIC telescopes is es- 
timated from Monte Carlo simulations. The simula- 
tion program of MAGIC is divided into three stages. 
The CORSIKA [3] program simulates the air showers 
initiated by either VHE 7-rays or hadrons. Here, we 
used the CORSIKA version 6.019, the EGS4 code for 
electromagnetic shower generation and QGSJET-II and 
FLUKA for high and low energy hadronic interactions 
respectively. The atmospheric model is based on studies 
of total mass density as a function of the height. The sec- 
ond stage of the simulation accounts for the Cherenkov 
light absorption and scattering in the atmosphere and 
then performs the reflection of these photons on the 
mirror dish. Finally, the last stage simulates the behavior 
of the MAGIC photomultipliers, the trigger system and 
the data acquisition electronics. Pulse shapes, noise 
levels and gain fluctuations obtained from the MAGIC-I 
data have been implemented in the MAGIC-I simulation. 
For MAGIC-II, results from laboratory measurement are 
implemented. 

For the study presented here, we simulated 3 x 10 6 7- 
rays with an energy distribution ranging from lOGeV to 
30 TeV and following a power law with a spectral index 
of -1.6, and impact point within 300 m to the center of 
the array. The simulated 7-ray source is point-like and 
located at the center of the camera. For the background 
simulation, we generated 2 x 10 7 proton showers (re- 
used 160 times) with an energy distribution ranging from 
30GeV to 30 TeV with a spectral index of -1.78 over a 
1200 x 1200m 2 area and a 5°-radius view cone. For 
both 7-rays and protons, the telescope pointing ranges 
from 5° to 30° in zenith distance and from 0° to 360° 
in azimuth. 

The energy distributions have been chosen harder 
than the typical 7-ray sources in order to increase the 
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Fig. 1 : Angular resolution of the MAGIC observatory in 
single mode (MAGIC-I) and in stereoscopic mode. 
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Fig. 2: Mean error on the impact parameter as a function 
of the impact parameter for different size cuts. 



number of event at high energy. Nevertheless, in the 
analysis, each event is weighted as a function of its 
energy to reproduce the spectral index of a realistic 
source (~— 2.5 for the Crab nebula) and of the CR back- 
ground (~— 2.7 for protons). The analysis of simulated 
events is performed with standard MAGIC Analysis and 
Reconstruction Software (MARS) Jffl which can analyze 
both single-telescope and stereo data. 

III. Analysis and Results 

A. Image parametrization 

Cherenkov images of each telescope are calibrated, 
extracted and parameterized independently. The algo- 
rithms used here, are those developed for the MAGIC-I 
data analysis. The image cleaning method is the two- 
threshold method using time information [5 1 with a high 
threshold of 6 photo-electrons within a time window of 
±4.5 ns, and a low threshold of 3 photo-electrons with 
time constraint of 1.5 ns. For each image, a set of param- 
eters (Hillas parameters |6| and others) is obtained. In a 
IACT array, the Hillas parameters of multiple images 
of the same shower are generally combined in order 
to reduce the number of parameters. However, here we 
preferred to keep the parameters separated as they are 
only twice and the 2 cameras have different designs. 

B. Angular resolution 

In single-telescope mode, the shower direction is 
reconstructed with the DISP method Q which deter- 
mines the distance between the image centroid and the 
shower direction as a function of the image parameters. 
Discrimination between the two possible directions at a 
distance DISP from the centroid along the major image 
axis is performed thanks to the image asymmetry. With 
two telescopes, the primary direction and the ground 
impact point can be reconstructed by stereoscopy with 
the major axis of the 2 images (intersection of the two 
planes represented by the axes of the two images). 

There are several ways to define the angular reso- 
lution of a 7-ray observatory. Here, we consider the 
radius of the circle centred on the simulated source, 



containing 68% of the reconstructed events. Only events 
reconstructed at less than 0.4° are considered in this 
calculation. The others are considered as lost in the 
background. With the DISP method most of these events 
correspond to the case in which the sense of prop- 
agation has been wrongly determined (^20%). Thus, 
the wrongly oriented events do not affect the angular 
resolution, although they result in a reduction of the 
effective collection area after cuts. In stereoscopic mode, 
only events with nearly parallel images can be lost. 

Figure 1 shows the 68 %-containment radius for 7-rays 
as a function of the energy for the DISP method with 
MAGIC-I data and for the stereo reconstruction method. 
The performance of the stereo reconstruction depends 
on the angle between the two shower-image axes. For 
angles below ^30°, the DISP method is generally better 
than the stereo reconstruction. Selection of events with 
an angle between the shower-image axes above 30° 
was done for the stereo angular resolution curve shown 
Figure 1 . Our analysis is still under development but the 
final reconstructed direction would be a combination of 
the DISP method applied on MAGIC-I and MAGIC - 
II data and the stereo reconstruction. Here we simply 
use the stereo reconstruction for events with an angle 
between image axes above 30° and the DISP method 
for the others. 

C. Energy resolution 

The total number of photo-electron in an image 
(image size) is strongly related to the primary energy 
and to the distance between the shower axis and the 
telescope (impact parameter). The energy reconstruction 
performance depends strongly on the impact parameter 
reconstruction. The stereoscopic method allows to re- 
construct the ground impact point (Figure 2) and also 
the altitude of the shower maximum (H max ) related to 
the distance between the shower direction and the image 
centroids. In stereo-data analysis, for each telescope an 
energy is reconstructed using lookup tables based on 
the image size, impact parameter, H max and distance 
to zenith. Combination of these two energy estimations 
provides the final reconstructed energy as well as a 
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(a) Case I 
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Fig. 3: Energy reconstruction of stereo data with lookup tables. Left-hand and right-hand panels show respectively 
the energy bias and the energy resolution as a function of the estimated energy. In the second plot, energy resolution 
obtained with stand alone MAGIC-I (RF method) is also shown for comparison. 



parameter describing the compatibility between these 
two energies. 

In single-telescope mode, the impact parameter and 
Umax cannot be properly reconstructed. The energy 
reconstruction is generally performed with the Random 
Forest (RF) technique [7| which is an optimized event 
classification tool trained with a subset of 7-ray simula- 
tion. 

Figure 3b compares the performance of the energy 
reconstruction of 7-rays with the lookup-table method 
applied on the stereo data and with the RF method using 
MAGIC-I data. Although the lookup-table method is 
simpler than the RF method, it provides a better energy 
resolution because the 3D-parameters of the shower are 
well recontracted by stereoscopy. The improvement is 
particularly interesting at low energies. A RF method 
using stereo data is under development and should 
improve the energy resolution even more. The results 
presented Figure 3b are obtained without any assumption 
on the source position (for both curves). It is possible to 
improve the energy resolution in stand-alone mode by 
assuming the source position (this is generally done for 
the spectral analysis of point-like source with MAGIC- 
I). In stereoscopic mode, such assumption provides 
negligible improvement as the impact parameter and 
Hmax are already well reconstructed (see figure 2). 
Spectral analysis performance of extended or poorly 
localized sources will strongly improve with the stereo- 
observation. 

D. Gamma/hadron separation 

The 7-ray/hadron separation is performed with the 
RF technique. A RF is trained to distinguish 7-ray and 
proton events with a set of Monte Carlo simulations. 
After training, any event can be classified by the RF, 
providing a parameter, called hadronness, distributed 
between and 1 . Values closer to mean that the event 
is more 7-ray-like and values closer to 1 mean that 
it is more hadron-like. The cosmic -ray background is 




Fig. 4: Hadronness distribution for 7-rays and protons 
with a reconstructed energy above 100 GeV. 

rejected mainly by a cut on the hadronness. Figure 4 
shows the hadronness distribution for protons and 7- 
rays. As expected, the proton distribution peaks at 1 
whereas the 7-ray distribution peaks at 0. The quality 
factor of the selection with this cut (Q — e 7 / \f[t P ) 
with e 7 and e p respectively the fraction of 7-rays and 
protons passing the cut) is Q — 3.1 at 100 GeV and 
Q = 6.5 at 500 GeV. The other hadrons are more easily 
rejected than protons, so most of the cosmic rays passing 
the selection cuts are protons. The background rate can 
be fairly estimated with only proton simulations. 

E. Sensitivity 

The improvement of the 7-ray/hadron separation and 
energy reconstruction from the single-telescope mode to 
two-telescope system improves strongly the sensitivity 
and spectral analysis performance of MAGIC. The angu- 
lar resolution improvement also reduces the background 
for point-like source improving the sensitivity for this 
kind of source. Figure 5 shows the minimal integral flux 
of a point-like source above a given energy for a 5- 
standard-deviation detection, with an excess of at least 
10 events representing more than 5% of the residual 
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Fig. 5: Sensitivity of the MAGIC telescopes in stereo- 
scopic mode compared to the Crab Nebula spectrum and 
to the performance achieved by HESS. 

background, in 50 h of observation with the MAGIC 
telescopes in stereoscopic mode, as a function of this 
energy. The current version of MAGIC-I reaches at the 
best a sensitivity of 1.6% of the Crab nebula flux. In 
stereoscopic mode, MAGIC should achieve a sensitivity 
better than 1%. This sensitivity is comparable to the 
performance of VERITAS and HESS above 
150 GeV. Between 50 GeV and 150 GeV, MAGIC have 
the best sensitivity of the current IACT. Above few TeV, 
the MAGIC performance suffers of the relatively small 
trigger region of the MAGIC-I camera. 



In few months, the MAGIC telescopes will start 
observation in stereoscopic mode. Monte Carlo simu- 
lations of MAGIC show that the angular resolution, the 
energy resolution, the 7-ray/hadron separation and the 
sensitivity will strongly improve. In 50 h, MAGIC should 
detect point-like sources with a flux lower than 1 % of the 
Crab nebula flux. This sensitivity is comparable to the 
performance of other current IACT but with an energy 
threshold about 2 times lower (^60 GeV). Moreover the 
performances shown here are obtained with a prelimi- 
nary stereo-data analysis that may probably be improved. 
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